Microlensing Effects in Atmospheres of Substars 



L.A. Berdina, A.A. Minakov 



Institute of RadioAstronomy, Kharkiv, Ukraine 
lberdina@mail.ru 

The purpose of the present work is the study of focusing properties of atmospheres of substars that is 
necessary for adequate interpreting of observational data and for solving the inverse problem consisting in 
\ recovery parameters of 'microlenses' (substars) and sources (quasars). Amplification factor for a quasar 
q . image as projected onto the field of microlenses-substars was computed for optical and radio wavelengths. 

(N : 

^ ■ Introduction 

It was discovered from observations of the Helix Nebula with the Hubble Space Telescope, that there exists 
O ' an extended cluster of objects near the central star, which can be referred as substars. Masses of substars 
] turned out to be much less than the mass of the Sun, which is patently insufficient for the thermonuclear 
reactions to start in their interiors. However, the substars are surrounded by the dense gaseous envelopes 
' (— i ' and have plasma coronas. In observations of distant sources (quasars) through the dense fields of substars, 
O ' refraction of rays in their atmospheres and gravitational fields should be taken into account. The refraction 
can result in the lens effect, which will cause distortions of parameters of the observed quasars. On the one 
O ; hand, careful analysis of these distortions is necessary to interpret observational data adequately. On the 
other hand, the observed clusters of substars are of interest for researchers as they can be used for monitoring 
of the well-known microlensing effect. 
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Model representation of the media surrounding a substar 

To investigate the lensing effect of the 'atmosphere' determination of the refractive index of the medium 
in the neighborhood of a substar is the most necessary. But it should be taken into account that the gas and 
plasma medium near the substar are situated in a substar's gravitational field. In this case, the multiplicative 
formula for the refractive index n{ r ) is valid: 

n{~r) = n g {~r)n med {~r), (1) 



where n g {r ) is the effective refractive index for the gravitational field, and n me( i{ r ) is the actual refractive 
index of the medium. Observations of distant sources through the fields of substars are possible only at small 
refraction angles from the Earth. Taken this fact into account the analysis of the lens effect was carried out 
in the paraxial optics approximation. 

Estimations show that the gravitational field near substar can be considered to be weak and approximately 
spherically symmetric. In this case the effective refractive index for the gravitational field is determined as 
n g (pr) ps 1 + ^f, r > R, where r g = 2< ^ 2 M is the gravitational radius of a substar, G is the gravitational 
constant, c is the velocity of light, and R is the substar radius. The gaseous envelope and corona near substar 
can also be considered approximately as spherically symmetric. Model representations show, that the density 
of particles in the atmosphere of substar decreases quickly as the distance from its surface increases. In this 
case the refractive index of the atmosphere close to substar can be presented as n mec [ ps 1 + P(r), P{r) < 1. 
Taken this into account the total refractive index ([TJ) can be represented as 

n(r) ps 1 + (3{r) + ^ (2) 
r 

Analysis of the chemical composition of atmospheres of substars reveal that hydrogen, helium and a dense 
plasma corona are the main components [6]. Considering this we imply (3(r)to be a sum of two components 
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- the gaseous (3 n (r) and plasma /3 e (r) ones: (3(r) 
0e(r) = -[ 



CO » UJp. 



Here u;„ 
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= (3 e (r) + (3 n (r). The plasma component is determined as: 
e, m e are electron charge and mass, N e is electron spatial 

density in the corona. Neutral gas component can be estimated using the theory of the dielectric polarization 
of gases. It is known that permittivity of gas e = n 2 (r) is connected with polarization a n and concentration 



of the gas atoms N n (r) as e n (r) = 1 + j^N n {r) [5]. Hence (3 n 



)N n {r). 



Analysis of focusing of radiation from an extended source in the ray optics 
approximation 

We consider a problem of propagation of radiation from a distant source through the medium with a 
specified refractive index. A method of the phase screen can be applied to solve the problem. Distortions 
of beams from a source occur in a small region near substar. Outside this region the beams practically do 
not differ from their rectilinear asymptotes therefore it can be considered that the rays reach a lens plane at 

z = propagating along the straight lines. Then they are refracted through an angle 6 ( r ) and propagate 
again along the straight lines (z > 0) [1]. Passing towards an observer the rays can intersect forming regions 
of the enhanced energy concentration. In this case the effect of the medium can be represented as the effect 
of a thin corrector, where the wave gains additional phase incursion kL{ p ) and p is the impact parameter, 
that is the shortest distance from nonperturbed light beam to the center of substar. In the small-angle 
approximation of the geometrical optics, the eikonal L(p) is determined by integration of n( r ) along the 

oo > 

undisturbed beam L(p) = J n(z, p )dz. The angle of deflection of a beam 9 (p ) at the screen is connected 

— oo 

with the eikonal by a simple relationship 6 (p) = VyL( p) [I]. 

The extended source S is situated at the distance z = —z s from substar. The point element of a source 

dp s determined by the vector p s radiates a field U s (p s )dp s . The position of the observer is determined by the 
vector p (offset of the observer from lens axis). In the small-angle approximation the point element generates 
the field dU p (jp) in the point P of the lens plane. According to the Huygens principle the total field in the 
point of observation formed by all surface elements of the source will look like the following: 



U p (z p ) 



k 
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exp(ik(z s + z p )) I dp I U s (p s ) exp ik 



(~p - pi) 2 p 



2z, 



dp s - (3) 



In observations we deal with casually radiating sources, therefore the field has a random character 
at the observer's plane. The average intensity in the point of observation is determined as: (I p (z p )) = 

(Up(zp)u;(z p )). 

For the further evaluations the source model should be specified as a totality of incoherently radiating 
point elements: (U s (p' s )U* (p")) = IsPsHp's ~Ps)- Here I s (pl) is the intensity distribution over the source 



surface specified as I s (p s ) 



exp 



Ps, 

(ft-*)' 

2R1 



, where P s is the offset of the centre of radiation, Iq is 



the intensity of the whole source surface and R s is the effective radius of the source. For the given source 
model after integrating over p s we obtain the following expression for intensity: 



k 2 Io 
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dp / dpx 



x exp 



k 2 R 2 o 
^zf P+tk 



— + — J PP ~ —P + L(p+ Ip] - L (p - \p 



(4) 

The characteristic scale for changes of L (p) is large comparing with the effective region dimensions p e f 
which provides the basic contribution to the integral over p. In this case the following simplification 
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can be introduced in Q: L (p + 2P) — L (p — « VL (p) p = 6 (p) p. As a result, after integrating over p, 
we obtain: 



-oo 



where F(~p) = A ^"p* + z6(~p ) — ^-P S J , z = /^f ° ■ The integral j5]) is a Laplace type integral. The main 
contribution to the integral is from the region p located near the points, where the exponent equals to 0: 

-P s = ~p +z9(^) (6) 

In statistics the formula J5j) characterizes the Huygens principle for intensity and the equation {6|) is 
known as the equation of lens or the aberration equation. It is more convenient to analyse the integral 

(JSJ) in angular coordinates, which better suits to ground-based observations: tp = j- (& current angle of 

observation), ipo = - ^f- (the angular size of the source), ib s = p ° (the angular coordinate of a maximum 

of the source radiation). Thus, the amplification factor for the extended source defined as the ratio of the 
intensity integrated over the angles at the point of observation to the intensity, which would be observed in 
the absence of focusing, is determined as 

oo 

« = m exp {-Wo) J exp r ~w 1 10 ) **■ (7) 

Here Iq(x) is a modified Bessel function, and <& (ip) = 1 + z J~i' . 

To estimate integral J7J we use the method of asymptotic estimation of Laplace type integrals [3j and 
thus obtain: 



2^ 1 _____ ( € \ T ( € 



1 (^) = 7^7, r 7 x zr exp - I, 



2d>'(^) Vo "V #o/ V# 2 

The physical sense of the root ^ received from the equality <&(ip)ip = is that it is a radius of a ring 
in a lens plane, which will be seen by an observer locating strictly on the focal semi-axis. Analysis of the 
dependence of q on the value of shift ip s shows, that q(ip s ) reaches a maximum in two cases: when the source 
is projected near the lens axis (tp s = 0) and near the caustic (ip a = tp sc )- For a central source the maximum 

amplification is determined as q = 2 &(^)ip • The secon d critical region of the angles ip is determined from 

a requirement ^ [<3? (ip) ip] = 0. Substituting a root of the equation ij) = ipcr in the equation of lens we 

obtain the angular coordinate of a caustic in the plane of the source positions ip s = ip sc . Calculation of 
the maximum amplification in projection of the source strictly at a caustic gives the following quantity: 

Numerical calculations 

As we recede from the substar surface, concentration of neutral particles and electrons decreases according 

to the power law: N e (r) = A e o (£) " , N n (r) = N n Q f-^) " where N e o >n o are the values of concentration 
at r = R. In the case of the selected model parameters of substar the refractive index of the medium and 
the deflection angle of a beam are equal to 

n(r) = l-Oe-(£) +a -'(f) +a f(7)> W 
a i s \ 2 (P\~ h e a n (p\-h n rp 
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where a e ,a n ,a g - coefficients of the gaseous, plasma components and of gravitational field: a e = ^,a n = 

a " EB "° , a g = and also Aq = — r ( fee -j 2 ^ , ao = — n r(h ) 2 " g amma -f unc ti° n ) • F° r numerical 
estimations it is possible to select following values of parameters of a substar [6], [2]: /i e ,n ~ 100 -j- 1000, 
iV n o ~ 10 18 cm~ 3 , A^eo ~ 10 s cm~ 3 , with mass M s ~ 0.01 M and radius R ~ O.3i? 

The analysis of (E| and estimations of the source brightness amplification allows to conclude the following 

• the field of gravity and the gaseous component deflect the light beams aside a substar, while the electron 
component scatters them. However, taking into account that the electron component decreases with 
increasing p most quickly, the atmosphere converges beams for sufficiently large p, or in other words, 
works as a positive lens. On the other hand, gaseous component acts at small distances from substar 
centre, since the concentration of the gas particles in the atmosphere decreases as the impact parameter 
increases. Thus, at large p the gravitational field provides the main contribution; 

• in optics a neutral atmosphere and the field of gravity of substar provide a positive lens effect. This 
results in forming a region of the enhanced amplification behind the lens (a focal semi-axis). For quasars 
with the angular sizes 10 -9 the magnification value is about 10. In the radio frequencies the plasma 
corona forms a conical caustic surface, at which the amplification exceeds 1 insignificantly. Nevertheless, 
the caustic will exhibit itself at the radio frequencies because of the presence of a caustic shadow p] 

Conclusions 

Atmospheres and fields of gravity of substars can create the lens effect, amplifying radiation of distant 
sources - quasars, observed through the dense fields of substars. The continuous monitoring of the light 
curves of quasars can provide important information about mass distribution of substars, structure of their 
atmospheres, and also about the fine spatial structure of the emitting regions of quasars. 
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